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Abstract. Recent developments in the theory of weak decays of heavy flavours 
are reviewed. Applications to exclusive semileptonic B decays, the semileptonic 
branching ratio and charm counting, beauty lifetimes, and hadronic B decays are 
discussed. 

1 Introduction and theoretical concepts 

To discuss even the most significant recent theoretical developments in heavy- 
flavour physics in a single talk is a difficult task. Fortunately, this field is 
blooming and will continue to be of great importance in view of several new 
experimental facilities (B factories) to start operating in the near future. Be- 
low, f will review the latest theoretical developments in this field and discuss 
the most important phenomcnological applications. They concern semilep- 
tonic B decays and the measurements of the CKM parameters \ V c t\ and \ V u b\, 
the semileptonic branching ratio and charm yield in inclusive B decays, the 
lifetimes of beauty mesons and baryons, and hadronic B decays, including 
the rare decays into two light mesons. I will start with an introduction to the 
main theoretical concepts used in the analysis of these processes. 

The properties of hadrons containing a single heavy quark Q are char- 
acterized by the large separation of two length scales: the Compton wave 
length 1 jmQ of the heavy quark is much smaller than the typical size 1 //Iqcd 
of hadronic bound states in QCD (see the left plot in Fig. 1). In the limit 
vhq — ► oo, the configuration of the light degrees of freedom in the hadron 
becomes independent of the spin and flavour of the heavy quark, fn that 
limit there is a global SU(2n/j) spin-flavour symmetry of the strong inter- 
actions, where rih is the number of heavy-quark flavours [1]. This symme- 
try helps in understanding the spectroscopy and decays of heavy hadrons 
from first principles. It does not allow us to solve QCD, but to parametrize 
the strong-interaction effects of heavy-quark systems by a minimal number 
of reduced matrix elements, thus giving rise to nontrivial relations between 
observables. In particular, all form factors for the weak B — ► D^lv transi- 
tions are proportional to a universal function £(u>), where w — vb • is 
the product of the meson velocities. At the zero-recoil point w = 1, corre- 
sponding to vb = «£>(•), this function is normalized to unity: £(1) = 1. The 
symmetry-breaking corrections to the heavy-quark limit can be organized in 
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an expansion in powers of the small parameters a s (mQ) and ^qcd/tiq- A 
convenient way to do this is provided by the heavy-quark effective theory 
(HQET), whose purpose is to separate the short- and long-distance physics 
associated with the two length scales, making all dependence on the large 
mass scale mg explicit. This allows us to derive scaling laws relating differ- 
ent observables to each other. The philosophy behind the HQET is illustrated 
in the right plot in Fig. 1. 
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Fig. 1. Length scales of a heavy hadron, and construction of the HQET 



The effective Lagrangian of the HQET is 
£eff = h iv-Dh 



2itiq 



where v is the velocity of the hadron containing the heavy quark, and h is the 
heavy-quark field. The Lagrangian incorporates the spin-flavour symmetry 
to leading order in 1/toq. The only nontrivial short-distance coefficient has 
recently been calculated to next-to-leading and even next-to-next-to-leading 
order [2] -[4]. The result is 

C(m Q ) = [a s (m Q )] 9/25 {l + 0.672a s (m Q ) + (f .33 ± 0M)a 2 s (m Q ) +...}. 

The long-distance physics is encoded in the hadronic matrix elements of 
HQET operators, e.g. 



£{H) = -{H\h{iDfh\H) , &(H) = 



(H\ha llv G'"'h\H) 



A{H) = M H - m c 



2M H 



which are referred to as the "kinetic energy", the "chromomagnetic interac- 
tion" , and the "binding energy" , respectively. They are important parameters 
entering many applications of the heavy-quark expansion, the calculation of 
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inclusive decay rates and spectra in particular. Certain combinations of these 
parameters can be extracted from spectroscopy: 

H 2 G {B) = | (Ml, - Ml) « 0.36 GeV 2 , n 2 G (A b ) = , 



nl(A b )-nl(B) 



M B M D f M Ab -M Aa _ 3 Mb* - M D - _ 1 
M B - Mb 4 M B - M D 4 



A(A b )-A(B) w 0.31 GeV. 

The individual parameters /zj and /l are scheme dependent, however, because 
of renormalon ambiguities in their definition. Once a scheme is chosen, they 
can be calculated nonperturbatively, or extracted from moments of inclusive 
decay spectra. Table 1 shows a collection of some experimental determina- 
tions of these parameters in the on-shcll scheme. 



Table 1. Determinations of the parameters A and fj% from inclusive decay spectra 



Reference 


Method 


A{B) [GeV] fil{B) [GeV 2 ] 


Falk ct al. [5] 


Hadron Spectrum 


fa0.45 « 0.1 


Gremm et al. [6] 
Chernyak [7] 
Gremm, Stewart [8] 


Lepton Spectrum 
(B -> Xlv) 


0.39 ±0.11 0.19 ±0.10 
0.28 ±0.04 0.14 ±0.03 
0.33 ±0.11 0.17 ±0.10 


Li, Yu [9] 


Photon Spectrum 
(B -> X Bl ) 


U - DO -0.30 u -' 1 -0.70 



Using the operator product expansion, any inclusive decay rate of a beauty 
hadron can be expanded as 

where {O n )n are the matrix elements of local four-quark operators, which 
parametrize nonspectator effects in these decays, and Cj are calculable short- 
distance coefficients, which depend on CKM parameters, the ratios of quark 
masses, and the renormalization scheme. The free quark decay emerges as the 
leading term in a systematic l/mj expansion, with bound-state corrections 
suppressed by two inverse powers of the heavy-quark mass. Note that ratios 
of inclusive decay rates are independent of /i 2 and the common factor mjj, as 
well as of most CKM parameters. The application of the operator product 
expansion to the calculation of inclusive decay rates relies on the assumption 
of quark-hadron duality. Strictly speaking, the theoretical description of such 
processes is thus not entirely from first principles. 
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2 Exclusive semileptonic decays 

The most important applications of the HQET concern the description of 
exclusive semileptonic decays based on the quark transition b — > civ. This is 
where the theory is well tested, and theoretical uncertainties are best under- 
stood. The most important result is a precision determination of the CKM 
parameter \V c b\ [10]. 

2.1 Determination of \V c b\ from B — ► D^lu decays 

The differential semileptonic decay rates as a function of the kinematical 
variable w = vb • v D ( t) are given by 

dr(B^D*iv) G 2 F M% N2 i— , 

- - F B -rl{l~r st ) 2 \Jw 2 - l(w+ l) 2 



dw 487T 3 

Aw 1 — 2wr* + rl 

1 + 



u> + l (1-r*) 2 

&r{B^Dlv) _ G 2 F M% 



\V cb \ z T\w) 



B „3 



dw 48ir 3 



r 3 (l + r) 2 {w 2 -lfl 2 \V cb \ 2 Q\w) 



where ru) = M D («) /Mb- In the heavy-quark limit, the form factors F{w) and 
Q(w) coincide with the universal function £(w) and are thus normalized to 
unity at w = 1. Much effort has gone into calculating the symmetry-breaking 
corrections to this limit, with the result that [11] 

= 1 + c A (a s ) + + S 1/m 2 + ... = 0.924 ± 0.027 , 

Q{\) = 1 + c v (a s ) + S' 1/m + S' 1/m2 + ... = 1.00 ± 0.07 , 

where the short-distance coefficients ca and cy are known to two-loop order 
[12, 13]. These numbers include the leading-logarithmic QED corrections. The 
absence of first-order power corrections to is a consequence of Luke's 

theorem [14]. The theoretical errors quoted above include the perturbative 
uncertainty and the uncertainty in the calculation (and truncation) of power 
corrections, added in quadrature. If instead the errors are added linearly, 
the result for .F(l) changes to 0.924 ± 0.041. This value has recently been 
confirmed in a different regularization scheme, in which the separation be- 
tween short- and long-distance contributions is achieved by means of a hard 
momentum cutoff [15]. 

A value for |V c f>| can be obtained by extrapolating experimental data for 
the differential decay rates to the zero-recoil point, using theoretical con- 
straints on the shape of the form factors. Model-independent bounds on the 
physical B — ► D^iD form factors can be derived using analyticity proper- 
ties of QCD correlators, unitarity and dispersion relations [16, 17]. Combining 
these methods with the approximate heavy-quark symmetry, very powerful 
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one-parameter functions can be derived, which approximate the physical form 
factors in the semileptonic region with an accuracy of better than 2% [18]. 
For instance, the function Q(w) can be parametrized as 

« 1 - Zp\z(w) + (51.pl - 10.)z 2 (w) - (252.pl - 84.)z 3 (w) , 

where z(w) = (\fw + 1 — \/2)/(^w + 1 + \/2), and p\ is the (negative) slope 
of the form factor at zero recoil. A similar parametrization can be given for 
the function T(w). At present, these constraints are only partially included 
in the analyses of experimental data. 

The world-average results of such analyses are [19, 20] 

\V cb \ T(l) = (35.2 ± 2.6) x 1CT 3 , 
1^10(1) - (38.6 ±4.1) x nr 3 . 

A good fraction of the present errors reflect the uncertainty in the extrapo- 
lation to zero recoil, which could be avoided by implementing the dispersive 
constraints mentioned above. When combined with the theoretical predic- 
tions for the normalization of the form factors at zero recoil, the data yield 
the accurate value 

|V c6 | = (38.2 ± 2.3 oxp ± 1.2th) x 1(T 3 , 

which is in good agreement with an independent determination from inclusive 
B decays [1, 19]. 



2.2 Tests of heavy-quark symmetry 

In general, the decays B D iv and B — > D*£D are described by four 
independent form factors: Q(w) for the former process, and h,Ai{w), Ri(w), 
B-2(w) - a combination of which defines the function T(w) for the latter 
one. In the heavy-quark limit, G(w), T(w) and h.Ai{w) become equal to the 
function £(w), whereas Ri(w) and R2(w) approach unity. The universality of 
the function £(w) can be tested by measuring the ratio Q(w) / T(w) as a func- 
tion of w. The ALEPH data for this ratio [21] are shown in Fig. 2; a similar 
measurement has also been reported by CLEO [22]. Within errors, the data 
are compatible with a universal form factor. At large recoil, where the exper- 
imental errors are smallest, this provides a test of heavy-quark symmetry at 
the level of 10-15%. 

A more refined analysis of symmetry-breaking corrections has been done 
by measuring the ratios Ri and i?2 close to zero recoil. There the HQET 
predicts that [1] 
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Fig. 2. Ratio of the two form factors Q(w) and T{w), from Ref. [21] 



4a s (m r ) /l 
J2i w 1 + 1 + - — = 1.3 ± 0.1 



3tt 



2m c 



i?2 



2m c 



0.8 ±0.2, 



/4 - = + — ^ = 0.2 ± 0.1 . 



G 



3(1 -n) 



The CLEO data for these three quantities are Ri = 1.24 ±0.29, R 2 = 0.72 ± 
0.19, and p A — p% = 0.20±0.19 [23], in good agreement with the theoretical 
predictions. With a little more precision, the data would start to test the 
pattern of symmetry-breaking effects in the heavy-quark expansion. 



2.3 Total exclusive semileptonic rates 

Summing up the average experimental results B(B — ► Dlv) = (1.95±0.27)%, 
B(B -» D*lv) = (5.05 ± 0.25)%, B(B D^hiD) = (2.30 ± 0.44)%, and 
adding B(B — > X u lv) = (0.15 ± 0.10)% for the inclusive branching ra- 
tio for semileptonic decays into charmless final states, one gets a total of 
(9.45 ±0.58)% [19], which is not far below the average value for the inclusive 
semileptonic branching ratio, B(B — > XI v) — (10.19 ± 0.37)%, measured at 
the T(As) resonance (see below). Thus, there is little room for extra contri- 
butions. 

A solid theoretical understanding of B decays into p-wave charm meson 
resonances would be important in order to address the question of whether 
there are additional contributions not accounted for in the above sum, and 
also to understand the main source of background in the determination of 
| V c b\ from B — ► D*£ v decays. The description of these processes in the context 
of the HQET involves, at leading order, two new universal functions: r 3 / 2 (w) 
for the decays into the narrow states (D\, D^), and 71/2(10) for those into the 
broad states (Dp, D\). In the heavy-quark limit, the differential (in w) decay 
rates vanish at zero recoil, since in that limit the p-wave states are orthogonal 
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to the ground state. There is a sizable l/m c correction to the B — > D\iu 
decay rate at zero recoil, which can be calculated in a model-independent 
way in terms of known charm meson masses [24]. It is important, since the 
kinematical region is restricted close to zero recoil (1 < w < 1.3). Detailed 
theoretical predictions for the scmilcptonic and hadronic decay rates into 
p-wave charm states, which incorporate the constraints imposed by heavy- 
quark symmetry, can be found in Refs. [24, 25]. 

3 Semileptonic b — > u decays and \V u b\ 

Two exclusive semileptonic B decay modes into charmless hadrons have been 
observed by CLEO; the corresponding branching ratios are [26] 

B(B^tt£v) = (1.8 ± 0.5 exp ± 0.2 model ) x Hr 4 , 

B(B -+plv) = (2.5i°l xp ± 0.5 modc i) x lCT 4 . 

That the theoretical description of these processes involves heavy-to-light 
form factors implies a certain amount of model dependence, since heavy- 
quark symmetry does not help to fix their normalization in a precise way. To 
some extent, a discrimination between models can be obtained by requiring 
a simultaneous fit of both exclusive channels. From a x 2 -weighted average of 
models, CLEO obtains [26] 

\V uh \ = (3.3 ± 0.4 cxp ± 0.7 modc i) x 1(T 3 . 

Even with the present, very limited statistics of the measurements, the the- 
oretical uncertainties are the limiting factor in the determination of \V u b\. In 
the future, the model dependence can and will be reduced mainly by combin- 
ing different theoretical approaches, in particular: lattice calculations, which 
are restricted to the region of large q 2 [27]; light-cone QCD sum rules, in- 
cluding 0(a s ) and higher- twist corrections [28, 29]; analyticity and unitarity 
constraints [30, 31]; dispersion relations [32]. Although the optimal strategy 
is not yet clear at present, I believe it will be possible to reach the level of 
15% theoretical uncertainty. 

The traditional way to extract |V^,| from inclusive B — > X Iv decays 
has been to look at the endpoint region of the charged-lepton energy spec- 
trum, where there is a tiny window not accessible to B — > X c I v decays. 
However, this method involves a large extrapolation and is plagued by un- 
controlled (and often underestimated) theoretical uncertainties. A reanalysis 
of the available experimental data (using the ISGW2 model) gives \V u b\ = 
(3.7±0.6 oxp ) x 10~ 3 [33], in agreement with the value quoted above. A better 
discrimination between b — > u and b — > c transitions should use vertex infor- 
mation combined with a cut on the invariant mass Mh (or energy Eh) of the 
hadronic final state [34]-[40]. Parton model calculations (with Fermi motion 
included) indicate that about 90% of all B — > X U £D decays have Mh < Md, 
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as shown in Fig. 3. Ideally, this cut would thus provide for a very efficient 
discriminator. In practise, there will be some leakage so that presumably one 
will be forced to require < M max with some threshold M max < Md- The 
task for theorists is to calculate the fraction of events with hadronic mass 
M h < M max , 

M max 



To calculate this fraction requires an ansatz for the "shape function" , which 
describes the Fermi motion of the b quark inside the B meson [41, 42]. The 
first three moments of this function are determined in terms of known HQET 
matrix elements. Still, some theoretical uncertainty remains, mainly associ- 
ated with the values of the 6-quark mass and the kinetic energy as well 
as unknown 0{a 2 s ) corrections. As an example, Fig. 3 shows the dependence 
of <P(M max ) on the value of m b . It turns out that the resulting theoretical 
uncertainty strongly depends on the value of the threshold M max . First esti- 
mates yield 5\V ub \/\V ub \ w 10% for M max = M D , and S\V ub \/\V ub \ w 20% for 
-Mmax = 1.5 GeV [38, 39]. This new method is challenging both for theorists 
and for experimenters, but it is superior to the endpoint method. I believe 
that ultimately a theoretical accuracy of 10% can be reached. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

M h [GeV] M max [GeV] 



Fig. 3. Hadronic mass distribution (left) and fraction of B — > X u £i/ decays with 
invariant hadronic mass below M max (right), from Ref. [38] 



4 Semileptonic branching ratio and charm multiplicity 

For many years, the apparent discrepancy between (some) measurements of 
the semileptonic branching ratio of B mesons and theoretical predictions for 
this quantity has given rise to controverse dispute and speculations about 
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deviations from the Standard Model. Another important input to this dis- 
cussion has been the charm multiplicity n c , i.e. the average number of charm 
(or anticharm) quarks in the hadronic final state of a B decay. From the 
theoretical point of view, 

SL Ee r{B^xiD) + /^d + r raxe ' 

n c = 1 + B(B -> X cS ) - B(B -> no charm) 

are governed by the same partial inclusive decay rates. The theoretical predic- 
tions for these two quantities depend mainly on two parameters: the quark- 
mass ratio m c /mb and the renormalization scale [i [43, 44]. The latter depen- 
dence reflects our ignorance about higher-order QCD corrections to the decay 
rates. The results obtained by allowing reasonable ranges for these parame- 
ters are represented by the dark-shaded area in Fig. 4. The two data points 
show the average experimental results obtained from experiments operating 
at the r(4a) resonance: B SL = (10.19 ± 0.37)% and n c = 1.12 ± 0.05, and at 
the Z resonance: B SL = (11.12 ±0.20)% and n c = 1.20 ±0.07 [19, 45]. At this 
conference, it has been emphasized that a dedicated reanalysis of the LEP 
data for the semileptonic branching ratio is necessary, because some sources 
of systematic errors had previously been underestimated [45] . To account for 
this, I have doubled the corresponding error bar on Bsl in Fig. 4. If we ignore 
the LEP point for the moment, it appears that the theoretical predictions for 
both Bsl and n c lie significantly higher than the experimental results. 

1.4 
1.3 
c° 1.2 
1.1 
1 

8 9 10 11 12 13 14 

B SL (%) 

Fig. 4. Theory [44] versus experiment for the semileptonic branching ratio and 
charm multiplicity 

At this conference, several new ingredients have been presented that shed 
some light on this problem. First, there are some nontrivial tests of the the- 
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ory. Tabic 2 shows details of the theoretical calculation in the form of the 
dominant partial decay rates normalized to B(B — > X C £D), for m c /nib = 
0.29 ± 0.03, using the results of Rcfs. [43, 44, 46]. (Note that the values for 
r cud' ar e larger than the value 4.0 ± 0.4 often used in the literature.) It has 
been argued that the assumption of local quark-hadron duality, which under- 
lies the theoretical treatment of inclusive decay rates, may fail for the decays 
B — > X cEs i , because there is only little kinetic energy released to the final 
state. To test this hypothesis, one can eliminate the corresponding partial 
decay rate, in which case one obtains a linear relation between n c and Bsl 
[47]. The result is shown as the light band in Fig. 4. The fact that, within 
this band, the original prediction (i.e. the dark-shaded area) is closest to the 
data indicates that there is no problem with quark-hadron duality. Another 
important check is provided by an experimental determination of the ratio 
r C ud' using flavour-specific measurements of charm branching ratios [48] . The 
result is [49]-[51] 

B(B — > open c) — B(B — > open c) . , „ , „ „ 

rend- = — % - (2 + r CTP - tw) = 4.1 ± 0.7 , 

£>SL 

in good agreement with the theoretical predictions given in Table 2. The 
theory input r CTP — r uSs > = 0.19 ± 0.03 in this extraction has a very small 
uncertainty. In summary, it appears that the heavy-quark expansion works 
well for both of the hadronic decay rates. 



Table 2. Predictions for ratios of partial inclusive decay rates 







^ cud' 


faces' 


^no charm 


fj, — rah 


0.22 =f 0.03 


4.21 ±0.01 


1.89 =F 0.44 


0.14 ±0.04 


fj, — mb/2 


0.23 =F 0.03 


4.75 ±0.02 


2.20 =F 0.49 


0.19 ±0.04 



4.1 Is there a "missing charm puzzle"? 

Several suggestions have been made to explain why the experimental value 
of n c = 1.12 ± 0.05 measured by CLEO [19] is smaller than the theoretical 
prediction n c = 1.20 ± 0.06 [44]. I have discussed above that the problem 
cannot be blamed on violations of quark-hadron duality, as was originally 
speculated by some authors. An interesting new proposal made by Dunictz 
et al. [49] is based on the fact that the theoretical definition of n c refers to a 
fully inclusive quantity counting the number of charm and anticharm quarks 
per B decay, irrespective of whether they end up as "open" or "hidden" 
charm. If there were a sizable branching ratio for decays into hadronic final 
states containing undetected (cc) pairs, 



B(B — > (cc) undetected + X) = 6( cg ) , 



19: Heavy-Quark Effective Theory 



11 



then the experimentally observed value of n c would be lower than the theoret- 
ical one: n° hs = n* h — 2bt cS y Note that the conventional charmonium states 
J/i/), tp' , Xci, Xc2, Vc are included in the charm counting and thus do not 
contribute to 6( cg ). However, it has been speculated that there may exist an 
enhanced production of exotic (ccg) hybrids, which decay into noncharmed 
light mesons and could yield a value of bug) of about 10% [49, 52]. 

Alternatively, enhanced flavour-changing neutral current (FCNC) pro- 
cesses such as b — > sg could, simultaneously, lower the predictions for BgL 
and n c by a factor of [1 + B(new FCNC)] -1 , thus providing for a "new physics 
explanation" of the missing charm puzzle [53, 54]. 




Fig. 5. Mechanisms for D and D production in B-meson decays 



At least a partial answer to the question of whether there is any support 
for such speculations is provided by new, flavour-specific measurements of 
charm production. Using charm-particle-lepton correlations in tagged B de- 
cays, it is possible to measure the relative rate of the two mechanisms shown 
in Fig. 5. Whereas it was previously assumed that the first mechanism was re- 
sponsible for all Z?-meson production, recent measurements indicate a rather 
large "wrong-charm" yield [48]: 

B(B -»5l) 

-4- 7 = 0.100 ± 0.026 ± 0.016 , 

B(B —> D X) 

corresponding to a "wrong-charm" branching ratio B(B — > D X) = (7.9 ± 
2.2)%. This observation is supported by measurements of both exclusive [55] 
and inclusive [56, 57] production of DD meson pairs. This effect was not 
included in Monte Carlo -B-decay generators and may be partly responsible 
for the discrepancy in the CLEO and LEP measurements of Bsl and n c [49] . 
When combined with previously measured decay rates, the new CLEO result 
implies [48] 

B(B -> X cZ ) - b (c5) = (21.9 ± 3.7)% , 

which is close to the theoretical prediction B(B — ► X c£ ) = (22 ± 6)% [44], 
indicating that there is not much room for decays into undetected (cc) pairs, 
represented by 6( cg ) . This measurement is supported by a DELPHI result on 
the production of two open charm particles [58] . Furthermore, by measuring 
a double ratio of flavour-specific rates, CLEO obtains a bound on charmless 
modes that is largely independent of detection efficiencies and charm branch- 
ing ratios. The result is B(B -> no open charm) = (3.2 ± 4.0)% < 9.6% (90% 
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CL), which after subtraction of the known charmonium contributions implies 
[48] 

B(B -> no charm) + b (cS) = (0.2 ± 4.1)% < 6.8% (90% CL) . 

Given that the Standard Model prediction for the charmless rate is B(B — > 
no charm) = (1.6 ± 0.8)% [46], it appears that there is little room for cither 
hidden (cc) production or new physics contributions. Again, this conclusion 
is supported by a DELPHI result obtained using impact parameter measure- 
ments: B(B -> no open charm) = (4.5 ± 2.5)% < 8.4% (90% CL) [58]. 

Combining these two measurements, I conclude that 6( cg ) < 5% (90% CL) 
cannot be anomalously large. The same bound applies to other, nonstandard 
sources of charmless B decays, i.e. B(new FCNC) < 5% (90% CL). This con- 
clusion is supported by the DELPHI limit B(B — > X sg ) < 5% (95% CL) 
obtained from a study of the p± spectrum of charged kaons produced in B 
decays [58]. It must be noted, however, that a preliminary indication of a kaon 
excess at large p±, as expected from enhanced b — > sg transitions [50], has 
been reported by SLD at this conference [59] . Even though a definite conclu- 
sion can therefore not been drawn before these measurements become final, 
at present there is no compelling experimental evidence of any nonstandard 
physics in inclusive B decays. 

The two CLEO measurements quoted above can be combined to give a 
new determination of the charm multiplicity, in which the unknown quantity 
6( cg ) cancels out. The result 



1 



B(B — > X C c) — 6( C c) — B(B — > no charm) + br c $\ 



1.22 ±0.06 



is significantly higher than (though consistent with) the value 1.12 ± 0.05 
obtained using the conventional method of charm counting, and in excellent 
agreement with the theoretical prediction. I believe that, ultimately, this 
new way of measuring n c will be less affected by systematic uncertainties 
than the traditional one. It seem that the "missing charm puzzle" is about 
to disappear. 



5 Beauty lifetime ratios 

The current world-average experimental results for the lifetime ratios of dif- 
ferent beauty hadrons are [60] 

452 = 1-06 ±0.04, 4^4 -0.98 ±0.07, 4S = 0-78 ± 0.04 . 
t(B°) r(B d ) t(B°) 

Theory predicts that \t(B s )/t(B ( i) — 1| < 1% [61], and it will be very difficult 
to push the experimental accuracy to a level where one would become sen- 
sitive to sub-1% effects. The theoretical predictions for the other two ratios 
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have been analysed to third order in the heavy-quark expansion. Unfortu- 
nately, these predictions depend on some yet unknown "bag parameters" B\ , 
B 2 , £1, £2, B, r parametrizing the hadronic matrix elements of local four- 
quark operators. In terms of these parameters, the results are [44] 



2 , 

kiBt + k 2 B 2 + k 3 ei + k±e 2 



r 



= 0.98 + 16tt 2 L Bi + P2 B 2 + p 3 £i + Pi e 2 + (p 5 + p 6 B)r 

(±> u ) mi L J 



where fcj and Pi are short-distance coefficients, whose values depend on the 
ratio m c /mb and on the renormalization scale. The large- N c counting rules 
of QCD imply that Bi = O(l) and = 0(1/N C ). The factorization ap- 
proximation for the meson matrix elements suggests that Bi «J_ and £j « 
[62]. Similarly, the constituent quark model suggests that B m 1 and 
r w \ipg£ (0)| 2 /|'i/ , ^(0)| 2 - The parameter r is the most uncertain one enter- 
ing the predictions for the lifetime ratios. Existing theoretical estimates for 
this parameter range from 0.1 to 2. Some recent estimates can be found in 
Refs. [63, 64]. 

Without a reliable field-theoretical calculation of the bag parameters, no 
lifetime "predictions" can be obtained. However, one can see which ranges of 
the lifetime ratios can be covered using sensible values for the bag parameters 
[65]. To this end, I scan the following parameter space: Bi,B e [2/3,4/3], 
e l e [-1/3, 1/3], r e [0.25, 2.5], m c /m b = 0.29 ± 0.03, f B = (200 ± 20) MeV. 
The resulting distributions for the lifetime ratios are shown in Fig. 6. The 
most important observation of this exercise is that it is possible to repro- 
duce both lifetime ratios, within their experimental uncertainties, using the 
same set of input parameters. The theoretically allowed range for the ratio 
t(B~)/t(B°) is centered around the experimental value; however, the width 
of the allowed region is so large that no accurate "prediction" of this ratio 
could have been made. Any value between 0.8 and 1.3 could be easily ac- 
commodated by theory. The predictions for the ratio t(A] ] )/t{B°), on the 
other hand, center around a value of 0.95, and the spread of the results is 
much narrower. Only a small tail extends into the region preferred by exper- 
iment. Requiring that the theoretical results be inside the 2a ellipse around 
the central experimental values, I find that the value of the parameter E\ 
must be close to zero, £1 w —(0.1 ± 0.1), in agreement with the expectation 
based on factorization. In addition, this requirement maps out the region of 
parameter space for e 2 and r shown in the right lower plot of Fig. 6. Thus, 
the data indicate a large (with respect to most model predictions) value of r 
and a negative value of £2 . Ultimately, it will be important to perform reliable 
field-theoretical calculations of the hadronic parameters, for instance using 
lattice gauge theory or QCD sum rules, to see whether indeed the observed 
beauty lifetime ratios can be accounted for by the heavy-quark expansion 
(for a recent sum-rule estimate of the parameters Bi and Si, see Ref. [66]). 
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Fig. 6. Upper plots: Theoretical distributions for beauty lifetime ratios. The two 
curves in each figure correspond to different choices of the renormalization scale. 
The shaded areas show the experimental values with their 2a error bands. Lower 
left: Two-dimensional distribution of lifetime ratios, together with the la and 2a 
contours around the central experimental values. Lower right: Distribution of the 
hadronic parameters £2 and r for simulations with results inside the 2a ellipse. 



This is particularly important for the parameter r, large values of which are 
not excluded a priori. Using a variant of the quark-model relation derived in 
Ref. [63], combined with a preliminary DELPHI measurement of the mass- 
splitting between the and Sf, baryons [67] , I find 



4 M%. - Ml b 
3 Aff , - Mf 



1.8 ±0.5. 



With such a large value, it is possible to explain the short Af, lifetime without 
invoking violations of local quark-hadron duality. 

Although at present there is thus no convincing evidence that the low 
value of the Ab lifetime could not be accommodated in the context of the 
heavy-quark expansion, it remains a possibility that in the future such a dis- 
crepancy may emerge, for instance if lattice calculations would show that 
r < 1. Then, one would have to blame violations of local quark-hadron 
duality to be responsible for the failure of the heavy-quark expansion for 
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inclusive, nonleptonic B-decay rates. Recently, several authors have studied 
in QCD-inspircd models how quark-hadron duality may be violated. Based 
on a simple model for the difference of two heavy-light current correlators 
in the chiral and large-./V c limits, Shifman has argued that deviations from 
local duality are due to the asymptotic nature of the operator product ex- 
pansion (OPE) [68]. More recently, Blok et al. have studied QCD 2 (i.e. QCD 
in 1 + 1 space-time dimensions) at next-to- leading order in 1/N C , where finite 
resonance widths provide for a dynamical smearing of correlation functions 
[69]. They find that, in the Euclidean region, deviations from local duality 
are indeed due to the divergence of the OPE, and are exponentially sup- 
pressed with Q 2 . In the physical (Minkowskian) region, on the other hand, 
they find that two mass scales, A\ and A 2) are dynamically generated. For 
A 2 < q 2 < A 2 , there are strong oscillations of correlation functions. Global 
duality works, but local duality may be strongly violated. For q 2 > A 2 , the 
oscillations become exponentially damped, and local duality is restored. They 
conjecture that a similar pattern may hold for real, four-dimensional QCD. 
Similar results have also been obtained by Grinstein and Lcbed [70], and 
by Chibisov et al. studying the instanton vacuum model [71]. Quite differ- 
ent conclusions have been reached by Colangelo et al. , who have considered a 
more "realistic" version of Shifman's model, finding that global duality holds, 
but local duality is spoiled by 1/Q effects not present in the OPE [72] . This 
study gives some support for the conjecture by Altarelli et al. that there may 
be non-OPE terms of order A/tuq increasing the nonleptonic decay rates of 
heavy hadrons [73]. Such terms could easily explain the short A^ lifetime. 
However, I stress that, although the model of Ref. [72] is interesting, it must 
not be taken as an existence proof of non-OPE power corrections. Nature 
may be careful enough to avoid power-like deviations from local duality. 

6 Hadronic B decays 

The strong interaction effects in hadronic decays are much more difficult 
to understand, even at a qualitative level, than in leptonic or semileptonic 
decays. The problem is that multiple gluon exchanges can redistribute the 
quarks in the final state of a hadronic decay. As a consequence, some phe- 
nomenological assumptions are unavoidable is trying to understand these 
processes. 

6.1 Energetic two-body decays 

It has been argued that energetic two-body decays with a large energy release 
are easier to understand because of the colour-transparency phenomenon: a 
pair of fast-moving quarks in a colour-singlet state, which is produced in 
a local interaction, effectively decouples from long- wavelength gluons [74]. 
This intuitive argument suggests that nonleptonic amplitudes factorize to a 
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good approximation. The main strong-interaction effects are then of a short- 
distance nature and simply renormalize the operators in the effective weak 
Hamiltonian [75]. For the case of b — + cud transitions, for instance, 

H cS = ^ V cb V; d [aifi) (du)(cb) + C2(jm) (cu)(db)} + ..., 

where (du) — ^7^(1 — 75)^ etc. are left-handed, colour-singlet quark cur- 
rents, and Ci(m;,) w 1.1 and 02(1111,) w —0.3 are Wilson coefficients taking 
into account the short-distance corrections arising from the exchange of hard 
gluons. The effects of soft gluons remain in the hadronic matrix elements of 
the local four-quark operators. In general, a reliable field-theoretical calcu- 
lation of these matrix elements is the obstacle to a quantitative theory of 
hadronic weak decays. 

Using Fierz identities, the four-quark operators in the effective Hamil- 
tonian may be rewritten in various forms. It is particularly convenient to 
rearrange them in such a way that the flavour quantum numbers of one of 
the quark currents match those of one of the hadrons in the final state of 
the considered decay process. As an example, consider the decays B — ► Dir. 
Omitting common factors, the various amplitudes can be written as 

Abo^d**- = (c 1 + ^-)(D+n-\(du)(cb)\B°)+2c 2 (D+n-\(dt a u)(ct a b)\B°) 

= 01 (7r-|(du)|0) (D+\(cb)\B°) , 
A So ^ D o v o = (ca + (D°7r° I (cu) (db) \B°)+ 2c, (D"tt°\ (ct a u) (dt a b) \ B°) 

= a 2 (£>°|(cu)|0) (A(db)\B°) , 

where t a are the SU(3) colour matrices. The two classes of decays shown 
above are referred to as class- 1 and class-2, respectively. The factorized matrix 
elements in the last steps are known in terms of the meson decay constants 
/tt and fo, and the transition form factors for the decays B — > D and B — > n. 
Most of these quantities are accessible experimentally. Of course, the above 
matrix elements also contain other, nonfactorizable contributions. They are 
absorbed into the definition of the hadronic parameters a\ and a 2 , which in 
general are process dependent. Recently, some progress in the understanding 
of these parameters has been made, leading to the "generalized factorization 
hypothesis" that [76]-[78] 

«i«ci(mi,), a 2 ~ c 2 (m b ) + C,ci(m b ) , 

where Q is a process-independent hadronic parameter (for energetic two-body 
decays only!), which accounts for the dominant nonfactorizable contributions 
to the decay amplitudes. To derive these results, one combines the 1/N C 
expansion with an argument inspired by colour transparency [78]. 
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Table 3. Experimental tests of the generalized factorization hypothesis 







ai 


a 2 /ai 


B 


D { *\ 


1.08 ±0.06 


0.21 ±0.07 


B 




1.07 =F 0.07 


0.23 ±0.14 


B 






|a 2 | = 0.21 ±0.04 


B 




1.10 ±0.18 





Some experimental tests of the generalized factorization hypothesis arc 
shown in Table 3, where I quote the values of a\ and extracted from the 
analysis of different classes of decay modes, using data reported by CLEO 
[79]-[81]. Within the present experimental errors, there is indeed no evidence 
for any process dependence of the hadronic parameters. The generalized fac- 
torization prescription provides a simultaneous description of all measured 
Cabibbo-allowed two-body decays with a single parameter £ = 0.45 ± 0.05 
extracted from the data on class-2 decays [78]. So far, this theoretical frame- 
work is fully supported by the data. 

The factorization hypothesis can be used to obtain rather precise values 
for the decay constants of the D s and D* mesons, by comparing the theoret- 
ical predictions for various ratios of B° decay rates, 



r(D+ir-) ■ V U J ' r(D*+7T-) ■ V U 
r(D+ D ;-) = oj4 //d- \ 2 r(D*+p;-) = Lgg ( f D . 



r(D+ P -) • V f P J ' r(D*+ P -) ■ V f P 

with data. These predictions are rather clean for the following reasons: first, 
all decays involve class-1 transitions, so that deviations from factorization 
are probably very small; secondly, the parameter a\ cancels in the ratios; 
thirdly, the two processes in each ratio have a similar kinematics, so that the 
corresponding decay rates are sensitive to the same form factors, however 
evaluated at different q 2 values. Finally, some of the experimental systematic 
errors cancel in the ratios (however, I do not assume this in quoting errors 
below) . Combining these predictions with the average experimental branching 
ratios [79] yields 

f Ds = (234 ± 25) McV , f D . = (271 ± 33) MeV . 

The result for fo s is in good agreement with the value fo s = 250 ± 37 McV 
extracted from the leptonic decay D s — > p + v [82]. The ratio of decay con- 
stants, /d*//l> s = 1-16 ± 0.19, which cannot be determined from leptonic 
decays, is in good agreement with theoretical expectations [83, 84]. Finally, I 
note that, assuming SU(3)-breaking effects of order 10-20%, the established 



18 M. Neubert 



value of fjj s implies fn ^ 200 MeV, which is larger than most theoretical 
predictions. 

6.2 B Decays into two light mesons 

Since last year, CLEO has observed (or put strict upper limits on) a number of 
rare B decay modes, many of which have strongly suppressed tree amplitudes 
and are thus dominated by loop processes. The relevant quark diagrams for 
the transitions b — > s(d)qq are shown in Fig. 7, where I also give the powers of 
the Wolfenstein parameter A w 0.22 associated with these diagrams. Because 
the penguin digram shown on the right is a loop process, it is sensitive to 
new heavy particles and thus potentially probes physics beyond the Standard 
Model. 




Ms: 0(A 4 ) 0(A 2 ) for c,i 

b -4 d : 0(A 3 ) 0(A 3 ) 

Fig. 7. Tree and penguin diagrams for rare B decays into light hadrons 

Some of the experimental results that have caused a lot of excitement are 
(all in units of 10 -5 and averaged over CP-conjugate modes; upper limits at 
90% CL) [81, 85, 86]: 



B(B° - 


+ K+TT-) 


= 1.5±g;| ± 0.1 ± 0.1, 


B{B+ - 


-> K n+) 


= 2.3t^±0.3±0.2, 


B(B° - 


-> 7T + 7T~) 


= 0.7 ±0.4 < 1.5, 


B(B+ 


-> 7T+7T ) 


= 1.0±gj<2.0, 


B{B+ 


- K+rj') 


= 7.1±|i±0.9, 



B(B -> r)' + X s ) = 62 ± 16 ± 13 (2.0 < p n > [GeV] < 2.7) . 

In Table 4, 1 show the results of an analysis of these results in terms of SU(3)- 
invariant amplitudes representing diagrams having tree topology , colour- 
suppressed tree topology C^'\ penguin topology PW, and SU(3)-singlet pen- 
guin topology S^'\ where primed (unprimed) quantities refer to b — > s(d) 
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transitions [87]. The results of this analysis establish the presence of a sig- 
nificant SU(3)-singlet contribution S' active in the decays with an rj me- 
son in the final state. They also indicate a rather large "penguin pollution" 
\P/T\ w 0.3±0.1 in the decays B — > tttt. This has important (and optimistic) 
implications for CP-violation studies using these decay modes [88] . 

Table 4. Amplitude analysis for some rare B decays into two light mesons. The 
dominant amplitudes are highlighted. 



Decay mode 


SU(3) invariant amplitudes 


(in units of 10" 6 ) 


B° -> K+7T- 
B+ -> K°n+ 


-(T' + P>)) 
P> j 


\ 


\P'\ 2 a 16 ±4 
=> |P| 2 a 0.8 ±0.4 


B° -> Tr+TT" 
B+ -> TT+TT 


-(T + P) 1 
-^(T + C) j 


\ 


|T| 2 a 8 ±4 
=> |T' | 2 a 0.4 ±0.2 


B+ -> r,'K+ 


±(T' + C' + 3P' + 4S') 


|S"| 2 a 18 ±3 

no interference 

|S"| 2 a5±2 

constr. interference 



Several authors have discussed the question whether the penguin ampli- 
tudes in these transitions are anomalously large. There seems to be some 
agreement that a standard analysis, using next-to-leading order Wilson co- 
efficients combined with factorized matrix elements, can accommodate the 
experimental results for the exclusive rare decay rates [89, 90]. However, it 
has been stressed that significant long-distance contributions to the charm 
penguin are required to fit the data [91, 92]. This is not too surprising, since 
the energy release in B decays is such that the cc pair in the charm penguin 
is not far from its mass shell, and hence this penguin is really more a long- 
distance than a short-distance process [93]. This is illustrated in Fig. 8. In 
the case of B — > Kit decays, the box may represent final-state rescattering 
processes such as B — > D S D — > Kit [94]-[97]. If the final state contains an 
rj 1 meson, the box may represent an anomaly-mediated coupling of the rj' to 
glue, or an intrinsic charm component in the 77' wave function. At least to 
some extent these are different words for the same physics. 



. & 


strong 
QCD 
forces 







*- 


b s 


*- 





light hadrons 



Fig. 8. Long-distance effects in "charming penguins" 
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In a way, the charm penguin is all there is in b — > s FCNC processes [93] , 
because the unitarity of the CKM matrix implies that 

Kv qb p q « v*v cb (p c - p t ) . 

q—u,c,t 

The top penguin P t simply provides the GIM cutoff for large loop momenta, 
whereas the up penguin P u is strongly CKM suppressed. The construction 
of the effective weak Hamiltonian is used to separate the short- and long- 
distance contributions to the charm penguin, in a way that is illustrated in 
Fig. 9. The short-distance contributions from large loop momenta are con- 
tained in the local operators Q3,...,6 and Q$, while long-distance effects from 
small momenta are contained in the matrix elements of the current-current 
operators Qi,2- From the physical argument presented above, one expects 
sizable long-distance contributions, since the energy release in B decays is in 
the region of charm resonances. Therefore, purely short-distance estimates of 
the charm penguin may be very much misleading. 



c 




«=1,2 k < n (only sum is RG-invariant!) 

Fig. 9. Operator product expansion for the charm penguin 



6.3 Explanations of the r\' modes 

To explain the yield of rj mesons in rare B decays, in particular the large in- 
clusive 77' production rate, is a challenge to theorists. Currently, the situation 
is still controversial, and several proposals are being discussed. It has been 
suggested that the process b — > sg* followed by g* — ► rj g is enhanced, either 
by the anomalous coupling of the rj to glue [98]-[100], or because the initial 
b — > sg rate is enhanced by some new physics [101]. Because of its three- 
body nature, this mechanism would mainly be responsible for multi-particle 
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inclusive processes containing r{ mesons, but give a minor contribution to the 
two-body modes B — > r]'K^*\ It has also been argued that an enhanced rj' 
production could result from a significant intrinsic charm component in the 
rf wave function, either through the colour-singlet mechanism b — ► (cc) i s 
followed by (cc)i — > rj' [102, 103], or through the colour-octet mechanism 
b — > (cc)g s followed by (cc) 8 — > r/A [104], which is strongly favoured by the 
structure of the effective weak Hamiltonian. The relevance of the intrinsic 

(c) 

charm mechanism depends on the size of a "decay constant" f^, , estimates 
of which range from 6 MeV to 180 MeV. A clarification of this issue is needed 
(see, e.g., the discussion in Refs. [89, 105]). Yet other possibilities to explain 
the large rf yield have been explored in Refs. [106, 107]. 



Table 5. Model predictions for ratios of rare B decay rates 
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Cheng, Tseng 


Datta et al. 


Halperin et al. 
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0.36-0.38 


fa 2 


B+ 


-> n'K+ 




B+ 


-> r)K+ 


0.04-0.14 


0-0.02 


0.03-0.08 


0.02-0.09 




B+ 


T]'K+ 




B+ 




0.6-2.9 


0.5-8.2 


0.26-0.48 


1.0-1.4 




B+ 


V 'K*+ 





Future measurements of other decay channels, as well as of the distribu- 
tion of the invariant hadronic mass M(X S ) in the decay B — > r/'X s , will help 
to clarify the situation. Preliminary CLEO data indicate that relatively large 
values of M(X S ) are preferred [81, 108]. Another good discriminator between 
models is provided by the ratios of the various B — > rj^K^ decay rates. 
A collection of theoretical predictions for such ratios is shown in Table 5. 
In summary, the rare decays of B mesons into light mesons provide unique 
opportunities to extract many yet unknown hadronic matrix elements. At 
present, there are no convincing arguments that these processes could not 
be understood in a conventional Standard Model framework. In particular, 
because of the complexity of the dynamics of the rj' meson, I believe it is pre- 
mature to deduce evidence for new physics from the yield of rj' production in 
B decays. 

7 Summary and outlook 

Owing to the combined efforts of experimenters and theorists, there has been 
significant progress in heavy-quark physics in the past few years. The exclu- 
sive and inclusive semileptonic decays mediated by the transition b — > civ 
are understood from first principles, using the heavy-quark expansion. Heavy- 
quark effective theory (HQET) works well in describing these processes, and 
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starts being tested at the level of symmetry-breaking corrections. This, by 
itself, is a remarkable development, which puts the determination of \V c b\ = 
0.038 ± 0.003 on firm theoretical grounds. New ideas are being developed for 
getting a precise (~ 10%) determination of \V u f,\, combining exclusive and 
inclusive methods. The error in the current value, jV^I = (3.3 ±0.8) x 10~ 3 , 
is still dominated by the theoretical uncertainty. 

New flavour-specific measurements of semi-inclusive charm yields indicate 
that the "problem" of the charm deficit and low semileptonic branching ratio 
is disappearing. The theoretical predictions for B(B — > X c£ ) and B(B — > 
no charm) are confirmed by the data. There is no compelling evidence for 
violations of quark-hadron duality, nor for new physics or exotic decay modes 
in B decays. The low value of the At, lifetime remains a surprise, but can be 
accommodated in a small corner of parameter space. Lattice calculations of 
baryon matrix elements of four-quark operators would help to clarify the 
situation. However, this may also be the first sign of departures from local 
quark- hadron duality. 

Hadronic two-body decays of B mesons with a large energy release can be 
understood in terms of a generalized factorization hypothesis, which includes 
the leading nonfactorizable corrections. Finally, new data on B decays into 
two light mesons open a window to study the details of nonperturbative 
dynamics of QCD, including long-distance penguins and other exotica. So 
far, the Standard Model can account for the data, but we should be prepared 
to find surprises as measurements and theory become more precise. For now, 
it remains to stress that the observed penguin effects are large and thus raise 
hopes for having large CP asymmetries in these decays. Hence, there is more 
Beauty to come in a bright future. 
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